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CYCLOTRON  MASER  INSTABILITY  FOR  INTENSE  SOLID  ELECTRON  BEAMS 


Han  S.  Uhm 
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University  of  Maryland 
College  Park,  Maryland  20742 


and 

* 

Ronald  C.  Davidson 
Office  of  Fusion  Energy 
Department  of  Energy 
Washington,  D.  C.  20545 

ms 

The  cyclotron  maser  instability  for  a solid  relativistic 


electron  beam  propagating  parallel  to  a uniform  axial  magnetic 


field  is  investigated.  The  stability  analysis  is 

carried  out  within  the  framework  of  the  linearized  Vlasov- 
Maxwell  equations.  It  is  assumed  that  v/y<<l,  where  v is 
Budker's  parameter  and  ymc  is  the  electron  energy.  Stability 
properties  are  investigated  for  the  choice  of  equilibrium 
distribution  function  in  which  all  electrons  have  the  same 
value  of  total  perpendicular  energy,  the  same  value  of  axial 
velocity,  and  a step-function  distribution  in  canonical 
angular  momentum.  The  instability  growth  rate  is  calculated 
including  a determination  of  the  optimum  value  of  the  beam 
radius  Rq  for  maximum  growth.  It  is  found  that  the  maximum 
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growth  rate  for  a solid  beam  is  comparable  to  the  maximum 


growth  rate  for  a hollow  beam. 
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I.  INTRODUCTION 


1-9 

In  recent  years,  the  cyclotron  maser  instability  has  been 

extensively  investigated  with  particular  emphasis  on  the  Implications 

for  intense  microwave  generation.  For  the  most  part,  previous  theoretical 

analyses  of  this  instability  have  been  carried  out  for  hollow  electron 
6—9 

beams.  The  present  paper  examines  the  equilibrium  and  cyclotron 
maser  stability  properties  of  a solid  relativistic  electron  beam  within 
the  framework  of  the  linearized  Vlasov-Maxwell  equations,  including  a 
determination  of  the  optimum  value  of  the  beam  radius  Rq  for  maximum 
growth  rate. 

The  present  analysis  is  carried  out  for  an  infinitely  long 
electron  beam  propagating  parallel  to  a uniform  axial  magnetic  field 


Bofe 


Equilibrium  and  stability  properties  are  calculated  for  the 


specific  choice  of  electron  distribution  function  [Eq.  (3)], 


fe(H"P0’P2)  * 2^  «(K1-Y0«c2)«(P9-P0)«(Pi-^V2) 


where  Hx  is  the  perpendicular  energy,  is  the  canonical  angular 

momentum,  P is  the  axial  canonical  momentum,  6(x)  is  the  Heaviside 
z 

step  function,  and  n^,  y q,  y , Pq,  and  are  constants.  Equilibrium 
properties  are  examined  in  Sec.  II,  and  stability  properties  are 
investigated  in  Secs.  Ill  and  IV,  assuming  that  v/y<<l,  where  v is 
Budker's  parameter. 

Introducing  the  normalized  Doppler-shifted  eigenfrequency 

Z ■(m-kV-sui  )/id  , the  TE  mode  dispersion  relation  can  be  expressed 
s z c c 

as  [Eq.  (31)]. 


Zs  - ^ ,,/o0nJ2(<*0n 
Y 


>i2(Vr*Y28‘H.)'° 


3 


r 


1 


I 


< 


where  J^x)  is  the  Bessel  function  of  the  first  kind  of  order  £,, 

a>c“eBQ/Ymc  is  the  electron  cyclotron  frequency,  k Is  the  axial 

wavenumber,  u is  the  complex  eigenfrequency , is  the  nth 

2 2 2 

root  of  8j,"(Yq”1) /Yq»  s-1,2,3,...  denotes  the  magnetic 

harmonic  number,  and  Q and  H are  coupling  coefficients  [Eq.  (27)]. 

s s 

Equation  (31)  is  valid  when  the  group  velocity  of  the  vacuum  waveguide 
mode  is  approximately  equal  to  the  beam  velocity.  Evidently,  for 
given  s,  the  maximum  TE  mode  growth  rate  occurs  at  a value  of  Rq/Rc 

corresponding  to  Ro^Rc“a2l',otOn’  Here  a21  is  the  ^irst  root  of  J2^°21^"0" 

This  result  is  different  from  that  obtained  for  a hollow  electron 

. 8-9 

beam. 

A detailed  numerical  analysis  of  the  TE  mode  [Eq.  (21)]  and  TM  mode 

[Eq.  (22)]  dispersion  relations  is  presented  in  Sec.  IV  for  the  case 

where  the  electron  gyroradius  is  much  less  than  the  beam  radius,  l.e., 

r <<R-..  Two  features  are  noteworthy  from  the  numerical  analysis. 

L U 

First,  the  system  is  completely  stable  to  TM  mode  perturbations 
when  the  group  velocity  of  the  vacuum  waveguide  mode  is  equal  to  the 
bean  velocity.  This  appears  to  be  a property  unique  to  the 
equilibrium  distribution  function  in  Eq.  (3),  where  all  electrons  have 
the  same  axial  velocity.  Second,  for  TE  mode  perturbations,  the  maximum 
growth  rate  for  a solid  beam  is  comparable  to  that  of  a hollow  beam. 

In  this  context,  we  conclude  that  solid  relativistic  electron  beams 
may  also  be  effective  in  generating  intense  microwave  radiation  by  the 
cyclotron  maser  instability. 


r 


« 


I 


1 
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II.  EQUILIBRIUM  PROPERTIES  AND  BASIC  ASSUMPTIONS 

The  equilibrium  configuration  is  illustrated  in  Fig.  1.  It 
consists  of  an  unneutralized  electron  beam  that  is  infinite  in  axial 
extent  and  propagates  parallel  to  a uniform  axial  magnetic  field 
Bq^z.  The  radius  of  the  electron  beam  is  denoted  by  Rq,  and  a grounded 
cylindrical  conducting  wall  is  located  at  radius  r*Rc.  The  applied 
magnetic  field  provides  radial  confinement  of  the  electrons.  As 
shown  in  Fig.  1,  we  introduce  a cylindrical  polar  coordinate  system  (r,0,z). 
In  the  present  analysis,  we  assume  that  the  beam  radius  Rq  is  larger 
than  twice  the  thermal  electron  Larmor  radius  r^,  i.e.,  Rq>2tl.  It 
is  also  assumed  that 

v/y«l  (1) 

2 2 

where  v*Nee  /me  is  Budker's  parameter, 

Rc 

N »2n  [ dr  r n^(r)  (2) 

6 J0  6 

is  the  number  of  electrons  per  unit  axial  length,  n^(r)  is  the 

e 

equilibrium  electron  density,  c is  the  speed  of  light  in  vacuo,  -e 

* 2 

and  m are  the  electron  charge  and  rest  mass,  respectively,  and  ymc 
is  the  electron  energy  in  the  laboratory  frame.  The  inequality  in  Eq. 

(1)  indicates  that  the  beam  is  very  tenuous,  so  that  the  perturbed 
fields,  to  lowest  order,  can  be  approximated  by  the  vacuum  waveguide 
fields. ^ Consistent  with  the  low-density  assumption  in  Eq.  (1),  we 
also  neglect  the  influence  of  the  (weak)  equilibrium  self-electric 
and  self-magnetic  fields  that  are  produced  by  the  lack  of 
equilibrium  charge  and  current  neutralization. 


For  present  purposes,  we  assume  an  equilibrium  distribution  function 


5 


of  the  form. 


f“(H,Pe,Pz)  - ^ «(Hi-l'0»c2).(Pe-P0)S(rr^V2)  , (3, 

2 

where  n^-const.  is  the  electron  density  at  r**0,  Yq®c  is  the  electron 
energy  in  a frame  of  reference  moving  with  the  mean  axial  velocity 


Vz^z  of  the  electron  beam. 


P0~(./2c)<rK)B0 


is  the  minimum  canonical  angular  momentum  of  the  electrons,  and 

II  , x>0  , 

(5) 

0 , x<0  , 

is  the  Heaviside  step  function.  In  Eq.  (3), 


u 2 /u2  * 2 2.1/2  , 2 4.  2 2.1/2 
H-ymc  *(Ha+c  p ) -(m  c +c  £ ) 


is  the  total  electron  energy, 


pe"[rPe-(e/2c)rBo1 


is  the  canonical  angular  momentum,  and  P “p  is  the  axial  canonical 

z z 

momentum,  where  lower  case  £ denotes  mechanical  momentum  and  the 
equilibrium  self  fields  have  been  neglected  in  comparison  with  the 
applied  magnetic  field  Bq. 

Making  use  of  Eqs.  (3)  and  (4),  and  defining  the  effective 


electron  Larmor  radius  by 


v(v1)1/2c/v 


where  u •eBrt/mc  is  the  nonrelatlvistlc  electron  cyclotron  frequency, 
c 0 

it  is  straightforward  to  show  that  the  electron  density  can  be 
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expressed  as 


0,  V 

ne(r)*no 


0<r<ltO-rL  • 

2 „2 , 2 

r -R_+r, 


R0-rL<r<R0+rL 


Iisln-1QZL)  , 

2 n \ 2rrL  / ’ 

0 , otherwise  , 

which  is  illustrated  in  Fig.  2.  Substituting  Eq.  (9)  into  Eq.  (2) 
gives  the  number  of  electrons  per  unit  axial  length 


(9) 


Ne"irn0R0 


(10) 


Note  from  Eq.  (9)  that  the  parameter  introduced  in  Eq.  (4)  determines 
the  effective  beam  radius. 


I 
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III.  LINEARIZED  VLASOV-MAXWELL  EQUATIONS  FOR  A TENUOUS  BEAM 


In  this  section,  we  make  use  of  the  linearized  Vlasov-Maxwell 
equations  to  investigate  stability  properties  for  azimuthally  symmetric 
perturbations  (3/30-0)  about  a tenuous  electron  beam  equilibrium 
described  by  Eq.  (3).  We  adopt  a normal-mode  approach  in  which  all 
perturbations  are  assumed  to  vary  with  time  (t)  and  axial  coordinate 
(z)  according  to 


6i^(x,t)=ii)(r)exp[i(kz-a)t)  ] , 


where  lmw>0,  u>  is  the  complex  eigenfrequency , and  k is  the  axial 
wavenumber.  Consistent  with  Eq.  (1),  it  is  also  valid  in  lowest 

7 8 

order  to  approximate  the  perturbed  fields  by  the  vacuum  waveguide  fields. ’ 

Without  loss  of  generality,  we  assume  that  the  amplitudes  of  the  perturbed 

axial  electric  and  magnetic  fields  are  normalized  according  to  E (r-O)-l 

z 

and  Bz(r-0)-l.  After  some  straightforward  algebraic  manipulation,  we 
obtain 


8l,a0./Rcc 

‘SW’1 


IV 

r 

Jo 


dr  r Jl(o0nr/Rc)Je(r) 


(11) 


for  the  transverse  electric  (TE)  waveguide  modes,  and 


R 

^ 2 f C dr  r V60nr/R  ) 

[Vl(S0n)]  J° 


P Jr)  - JJr) 


kc 


(12) 


for  the  transverse  magnetic  (TM)  waveguide  modes.  [For  a detailed 
derivation  of  Eqs.  (11)  and  (12),  see  Ref.  8].  In  Eqs.  (11)  and  (12), 
J^(x)  is  the  Bessel  function  of  first  kind  of  order  £,  agn  and  Bgn 
are  the  nth  roots  of  J^(agn)-0  and  and  the  Pertur'>et^ 

charge  and  current  densities  are  defined  by 


8 


Pe(r)—  e|d3p  fe(r,?)  , (13) 

and 

£(r)=-ejd3p  ^ fe(r,^)  , (14) 

where  . The  perturbed  distribution  function  in  Eqs . (13)  and 

(14)  is  expressed  as 

r°  {-  3 o 

fe(r»S>“ej  dT  exp{i[(kpz/Ym)-a)]T}  <£(r’)  + > * g^r  , 

(15) 

where  T=t'-t,  and  the  particle  trajectories  * ( t * ) and  ^'(t')  satisfy 
d^'/dt'*^'  and  djj'/dt'^-e^’xBg^/c,  with  initial  conditions  £'(t'=t)»£ 
and  v'(t'“t)“v.  In  obtaining  Eq.  (15),  use  has  been  made  of  the  axial 


orbit 


z'-z+(pz/Ytn)  (t  '-t) 


For  present  purposes  it  is  assumed  that 


| ft  | “ | w-kVz-sa>  | <<u>c 


where  fis=ai-kVz-suc  is  the  Doppler-shifted  frequency,  u)c“eBg/Ymc 

is  the  electron  cyclotron  frequency  in  the  laboratory  frame,  and  s=l,2,3,. 

denotes  the  magnetic  harmonic  number.  To  simplify  the  present  analysis, 


we  also  assume  that 


a 2 
v/y  <<(BxwcRc/c)  , 


2 1/2 

where  Bx“(1-1/Yq)  • Equation  (17)  is  easily  satisfied  in  parameter 

regimes  of  experimental  interest. ^ ^ Within  the  context  of  Eq.  (17), 
it  is  valid  to  neglect  the  terms  proportional  to  3f^/3P0  in  Eq.  (15), 

2 

since  the  corrections  associated  with  these  terms  are  order  (v/ y) (c/BxwcRq) 
(<<1)  or  smaller. 


The  particle  trajectories  r'(T)  and  0'(t)  in  the  equilibrium  fields 
are  required  in  order  to  evaluate  the  perturbed  distribution  function 
in  Eq.  (15).  The  transverse  (radial  and  azimuthal)  motion  of  a typical 
electron  is  illustrated  in  Fig.  3(a)  for  a positive  value  of  canonical 
angular  momentum  (P  >0) , and  in  Fig.  3(b)  for  a negative  value  of 
canonical  angular  momentum  (P  <0) . (The  dotted  circle  is  the  electron 
orbit  in  a plane  perpendicular  to  the  z-axis) . The  radial  distance 
of  the  electron  from  the  z-axis  at  time  t'  = t and  at  t'=*t'  are  denoted 
by  r and  r' , respectively.  The  point  C is  the  gyrocenter  of  the 
trajectory.  The  angular  coordinates  <f>  and  <t>'  are  the  perpendicular 
velocity-space  polar  angles  at  times  t'^t  and  t'*t',  and  are 
related  by 

4>  ’ =*4>+(ujc/y)  (t  ’-t)  . 


The  transverse  velocities  at  time  t'=t  and  t’=t’  are  denoted  by  vx 
and  respectively,  and  the  corresponding  speeds  are  defined  by 
v1*(vr+vQ)  and  vj  = (V  +v^  ) 

The  time  integration  in  Eq.  (15)  can  be  carried  out  by  making  use 
of  the  Bessel  function  summation  theorem  for  the  triangles  OAB  and 
ABC  in  Fig.  3.  After  some  straightforward  but  tedious  algebra,  the 
perturbed  TE  mode  distribution  function  can  be  expressed  as 


eR  pA 


°0n 


cm 


/ kp  v 3f 


3f 
e 

3P. 


VW^cV 

yu-su)  -kp  /m 
c z 


l (i)S  exp[is'(0-4))Js,(«Onr/Rc)Js+s,(^^)  , 
I ■-«  c c 


(18) 


and  the  perturbed  TM  mode  distribution  function  is  given  by 
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f^(r^)-eY  — * 


3fe  Js(S0n^/n,b,cRc)  ? _s’+l 


3p  yw-sw  -kp  /m 


I U)! 


(19) 


»«pK.'(8-«JJ  .{B^r/R  )J  (^i)  , 

c c 

where  use  has  been  made  of  Vg»vxsin((t> '-9 ')  = (pA/ym)sin(<t> '-0 ')  and 
2 2 2 

px=Pr+pQ  [For  a detailed  derivation  of  Eqs.  (18)  and  (19),  the  steps  are 
similar  to  the  derivation  of  Eqs.  (38)  and  (39)  in  Ref.  8]. 

In  obtaining  Eqs.  (18)  and  (19) , we  have  neglected  the  mode  coupling 
between  different  values  of  s,  which  is  consistent  with  the  inequality 
in  Eq.  (16).  Contributions  from  the  perturbed  radial  electric  field 
Er(r)  have  also  been  neglected  in  Eq.  (19),  which  is  valid  provided 
the  group  velocity  of  the  vacuum  waveguide  mode  is  significantly 
different  from  the  beam  velocity.  [For  TM  mode  perturbations, 
however,  it  is  shown  in  Sec.  IV  that  the  system  is  stable  when 
the  group  velocity  of  the  waveguide  mode  is  approximately  equal  to 
the  beam  velocity. ] 

The  equilibrium  electron  distribution  function  in  Eq.  (3)  can  be 
expressed  as 


2-irmr  oi 
L c 


6 (pA-mr^  coc)  6 (P^-ymuj 


c>*  (! 


2 r.2  , 2 

r -R„+r, 


sin<(>  - 


~ 0 rL  \ 
2rr.  / 


(20) 


after  some  simple  algebraic  manipulation.  In  Eq.  (20),  $“$-0. 

Substituting  Eqs.  (18)-(20)  into  Eq.  (11)-(14),  and  making  use 
* 2 2 -1/2 

of  y^y^d-V^/c  ) and  Eq.  (10),  we  obtain  the  dispersion  relation 


'2  a > 

a)  , 2 On 

~2  -k  -JT 

\C  R / 

c 


4v 


f(»-kV2)Qs(a0nyRc,a0orL/Rc> 


*2 T , .2 
YRcJ2(a0n) 


to-kV  -so) 
z c 


(»2-k2c2)Hs  (»0nR0/Rc , vW  62' 


(21) 


(w-kV  -sw  ) ‘ 
z c 


for  TE  mode  perturbations,  and 
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*vV60.VVBO.rl/V 


u-kv 


ui-kV  -sin 
z c 


(kc-wV  /c)‘ 

7. 


(uj-kV  -sai  ) ‘ 
z c 


for  TM  mode  perturbations.  In  Eqs.  (21)  and  (22),  B1=r  u /c,  v* 

L C 

2 2 

Nge  /me  is  Budker's  parameter,  and  the  coupling  coefficients 

H , Q , and  K are  defined  by 
s s s 


(22) 


Hs(x,y)=(l-y/x)2J2(x-y)2J’(y)2+2J’(y)I^(x,y)  , 

Qs (x, y) -2 (1-y/x) 2J2 (x-y) 2 (y) 2+yJ^ (y) (y) ] 

(23) 

+2[2J’(y)+yJ"  (y) ]lf (x,y)+yJ’ (y) [I®  (x.y) 

3 S S S“1 

_Is+l(x,y) ^ ’ 

Kg  (x ,y)  - ( 1-y/x)  ^ (x-y)  2Jg  (y)  2+Jg  (y)  Ig  (x, y)  , 


E M 

where  the  prime  (')  denotes  d/dy.  The  integrals  I and  Ig  in  Eq.  (23) 
are  defined  by 


IE(X  y)= i — 

Vx’y;  . 2 

2irx 
xexp 


fx+y 


22  - 1/2 
du  u J,(u)j2  d4>sin<)>J  [(u  +y  -2uysin<j>)  ] 

0 S 


■ — J 

• 

j 

. -1 
iscos 

y-usini 

(u2+y2-2uysin^)1/2 

j 

tm,  , i ^ 

I„(x,y)  - — j 


irx 


<exp 


du  u Jq(u) 


« 2 2 » 1 / 2 
d<f>  Jg[(u  +y  -2uysin4>)  ] 


C a 

iscos 

y-usin<f 

, 2^  2 _ , -.1/2 
(u  +y  -2uysin<t>) 

(24) 


where 


. . 2^  2 2 . 

►o-81" 


For  the  case  of  small  thermal  electron  Larmor  radius, 


(25) 
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rL/R0«l,  (26) 

E M 

the  terms  proportional  to  1^  and  1^  in  Eq.  (23)  are  negligible  in 
comparison  with  other  terms.  We,  therefore,  approximate  Eq.  (23)  by 

Hs(x,y)=J2(x) 2jg(y)2  , 

Qs(x,y)-2J2(x)2[r (y)2+yj; (y)J^  (y) ] , (27) 

K (x,y)-J. (x)2J  (y)2  . 
s i s 


For  a specified  value  of  Budker's  parameter  (v/y)  and  perpendicular 

beam  energy,  it  is  evident  from  Eq.  (27)  that  the  coupling  coefficient 

between  the  vacuum  waveguide  mode  and  the  electron  cyclotron  resonance 

mode  (u=kV  +sui  ) is  a maximum  when  Jl(a„  Rrt/R  )=0  for  the  TE  mode,  and 
z c z Un  u c 

when  J,(8rt  R„/R  )=0  for  the  TM  mode.  In  this  context,  we  conclude 
1 Un  U c 

that  the  maximum  coupling  for  magnetic  harmonic  number  s occurs  for  a 
value  of  R_/R  given  by 

U C i 

a21/a0n  ’ TE  mode  • 


Rn/R  - 
0 c 


(28) 


“ll^On  ’ ™ mode  ’ 


where  and  are  the  first  roots  of  J2(ot2^)=0  and  J^(a^^)»0, 
respectively.  Note  from  Eq.  (28)  that  the  value  of  Rq/Rc  corresponding 
to  maximum  coupling  is  the  same  for  all  magnetic  harmonic  numbers, 

8 9 

which  is  remarkably  different  from  the  results  previously  obtained  ’ 

for  a hollow  electron  beam,  where  the  maximum  TE  mode  coupling  coefficient 

2 

is  proportional  to  Jg(as^)  • As  shown  in  Eq.  (27),  the  maximum  TE 

2 

mode  coupling  coefficient  for  a solid  beam  is  proportional  to  J2  (a^)  • 
Therefore,  for  purposes  of  producing  intense  microwave 
radiation,  it  appears  that  a solid  electron  beam  may  be  as  effective 
as  a hollow  beam.  In  the  following  section,  the  dispersion  relations 


"TT 


~ 
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in  Eqs.  (21)  and  (22)  are  solved  numerically  for  the  case  of  small 
thermal  Larmor  radius,  with  rL/RQ<<l. 


i 
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IV.  STABILITY  ANALYSIS 


In  this  section,  the  linear  dispersion  relations  in  Eqs.  (21) 
and  (22)  are  solved  numerically  for  rL<<RoI  Ror  t*'e  TE  mode  dispersion 
relation  [Eq.  (21)],  we  assume  that  the  group  velocity  of  the  vacuum 

g 

waveguide  mode  is  approximately  equal  to  the  beam  velocity,  i.e.. 


kc 

V - 

g u> 


Making  use  of  the  TE  vacuum  waveguide  mode  dispersion  relation, 

2 2 2 2 2 

w /c  -k  -a.  /R  , the  cyclotron  resonance  condition,  w=kV  +u>  , and  Eq. 
Un  c z c 

(29),  we  obtain  the  lowest-order  solution 


“O^z^c  ’ 


VYz8zS(VC  • 


R -a-  c/y  so)  , 
c On  z c 

2 2 2 -1 

for  microwave  generation  at  frequency  u)=y  sw  . Here  y *(1-8  ) 

z c z z 

and  Sz“vz/C-  Substituting  Eq.  (30)  into  Eq.  (21)  and  iterating,  we 
obtain  the  approximate  TE  mode  dispersion  relation 

2 

Zs  " "T1  L ))  tQS(x*y)VSYz^Hs(x*y)1“0  ’ (31) 

y On  2 On 

where  the  normalized  Doppler-shifted  eigenf requency  Zg  is  defined 

by  Z “(ai-kV-sa  )/ai  , Q and  H are  defined  in  Eq.  (27),  and  x-a„  Rrt/R 
' s z c c s s OnOc 

and  yman  rT /R  - y s8,.  Since  the  electron  beam  is  assumed  to  satisfy 
1 On  L c z i 

RQ>2rL  (see  Sec.  II),  the  allowable  range  of  rq/rc  is  restricted  to 

R_/R  >2y  s8,/a_  , where  use  has  been  made  of  Eq.  (29).  [Therefore, 
u c z * un 

in  Fig.  4 the  growth  rate  plots  are  presented  only  for  values  of 


VV2VB*/a0n-> 
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Shown  In  Fig.  4 are  plots  of  normalized  growth  rate 

versus  R-./R  obtained  from  Eq.  (31)  (TE  mode)  for  v“0.001,  B,“0.4, 
u c 

Bz“0.3,  and  several  values  of  (s,n).  Several  features  are  noteworthy 
from  Fig.  4.  First,  the  maximum  growth  rate  for  each  set  of  (s,n) 
occurs  at  the  value  of  rq/rc  given  by  Ro/^c“a2]/aOn*  Second,  the  general 
shape  of  the  growth  rate  curve  is  similar  in  form  to  that  obtained  for 
s“2  perturbations  about  a hollow  electron  beam.  [See  Fig.  7 of  Ref. 

8.]  Third,  the  maximum  growth  rate  for  the  solid  beam  configuration 
considered  here  is  comparable  to  the  maximum  growth  rate  for  a hollow 
beam.  **  In  this  context,  we  conclude  that  the  cyclotron  maser  instability 
in  a solid  beam  may  also  be  an  effective  means  for  generating  intense 
microwave  radiation. 


We  now  investigate  a typical  example  that  illustrates  the  dependence 

of  TM  mode  stability  behavior  on  k.  Figure  5 shows  a plot  of  normalized 

growth  rate  Imu/iD  “w./o>  versus  normalized  axial  wavenumber  kc/ui 
c i c c 

obtained  from  Eqs.  (22)  and  (27)  for  v-0.001,  R-/R  “0.5,  rT /R  “0.1, 

U C L C 


8 “0.3,  n“2  and  several  values  of  s.  It  is  evident  from  Fig.  5 and 
z 

2 

Eq.  (22)  that  the  growth  rate  vanishes  at  k-k^-y zBzSwc/c,  which  is 
identical  to  the  condition  in  Eq.  (29).  We  therefore  conclude 
that  the  system  is  completely  stable  to  TM  mode  perturbations  when 
the  group  velocity  of  the  vacuum  waveguide  mode  is  equal  to  the 
beam  velocity.  This  appears  to  be  a unique  property  of  the  equilibrium 
distribution  function  in  Eq.  (3),  where  all  electrons  have  the  same 
axial  velocity.  It  is  important  to  note,  however,  that  the  TM  mode 


growth  rate  can  be  substantial  for  axial  wavenumbers  different  from 
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V.  SUMMARY  AND  CONCLUSIONS 

In  this  paper,  we  have  examined  the  cyclotron  maser  instability 
for  an  intense  solid  electron  beam  propagating  parallel  to  a uniform 
axial  guide  field.  The  analysis  was  carried  out  within  the  framework 
of  the  linearized  Vlasov-Maxwell  equations.  The  equilibrium 
(Sec.  II)  and  cyclotron  maser  stability  (Secs.  Ill  and  IV) 
properties  were  investigated  in  detail  for  the  choice  of  distribution 
function  in  which  all  electrons  have  the  same  value  of  perpendicular 
energy,  the  same  value  of  axial  velocity,  and  a step-function  distribu- 
tion in  canonical  angular  momentum.  One  of  the  most  Important 
conclusions  of  this  study  is  that  the  maximum  instability  growth  rate 
for  a •’olid  electron  beam  is  comparable  to  that  of  a hollow  beam. 
Moreover,  it  is  also  found  that  the  maximum  growth  rate  for  magnetic 
harmonic  number  s occurs  at  a value  of  Rq/Rc  corresponding  to  Rg/Rc= 

a„  /an  for  TE  mode  perturbations,  and  Rn/R  for  TM  mode 

zi  un  u c un 

perturbations.  For  the  TM  mode,  however,  the  growth  rate  vanishes 
when  the  group  velocity  of  the  vacuum  waveguide  mode  is  equal  to  the 
beam  velocity.  In  this  context,  the  TE  mode  is  more  effective  at 
producing  intense  microwave  radiation  in  a solid  beam  than  the  TM  mode. 
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FIGURE  CAPTIONS 
— 

Fig.  1 Equilibrium  configuration  and  coordinate  system. 

Fig.  2 Electron  density  profile  [Eq.  (9)]. 

Fig.  3 Electron  orbit  for  (a)  P >0,  and  (b)  P <0  projected  on  a 
plane  perpendicular  to  the  z-axis. 

Fig.  4 Plots  of  normalized  TE  mode  growth  rate  versus  Rq/Rc 

[Eq.  (31)]  for  v=0.001,  8X=0.4,  8Z=0.3,  and  several  values 
of  (s,n). 

Fig.  5 Plots  of  normalized  TM  mode  growth  rate  ai. /w  versus  kc/u 

ice 

[Eq.  (22)]  for  v=0.001,  n=2,  8 =0.3,  R_/R  =0.5,  rT /R  =0.1, 

Z U C L C 

and  several  values  of  s. 
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